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ABSTRACT 

We study eigenmodes of acoustic oscillations of high multipolarity I ~ 100 - 1000 and 
high frequency (~ 100 kHz), localized in neutron star envelopes. We show that the 
oscillation problem is self-similar. Once the oscillation spectrum is calculated for a 
given equation of state (EOS) in the envelope and given stellar mass M and radius R, 
it can be rescaled to a star with any M and R (but the same EOS in the envelope). 
For I > 300 the modes can be subdivided into the outer and inner ones. The outer 
modes are mainly localized in the outer envelope. The inner modes are mostly localized 
near the neutron drip point, being associated with the softening of the EOS after the 
neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and 
accreted matter and show that the spectra of the inner modes are essentially different. 
A detection and identification of high-frequency pressure modes would allow one to 
infer M and R and determine also the EOS in the envelope (accreted or ground-state) 
providing thus a new and simple method to explore the main parameters and internal 
structure of neutron stars. 
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1 INTRODUCTION 

Neutron stars can be considered as resonators where vari¬ 
ous oscillation modes can be excited. These oscillations are 
attracting much attention because, in principle, they can 
be used to study the internal structure of neutron stars. 
Some of them (for instance, r-modes) can be accompanied by 
gravitational radiation. Because neutron stars are relativis¬ 
tic objects, their oscillations must be studied in the frame¬ 
work of General Relativity. The relativistic theory of oscil¬ 
lations was developed in a series of paper s by Thorne and 
coaut h ors ||TIloriie_J^Campolihtero| |l96j; |Price_fc_Thomj 
^G^Jrh^ Icampolattoro^^ThomX * ^^^ 
IlDser fc Thornell97.3ll . In particular, the rapid (~ 1 s) damp¬ 
ing of p-modes with multipolar ity I = 2 by g ravitational 
radiation was demonstrated bv iThorn^ dlQbQall . An exact 
treatment of general-relativistic effects is complicated, but 
in many cases it is possible to use the relativistic Cowling 
approximation iMcDermott et a illl98.1l . An analysis of var¬ 
ious oscillation mo des and mechanism s for their dissipation 
was carried out lw^^ennot^_^^8^ . Let us also note the 
review paper bv IStergioulaj (l200.3ll . which contains an ex¬ 
tensive bibliography. As a rule, one considers neutron star 
oscillations with low values of 1. 

Although neutron stars are objects at the final stage of 
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stellar evolution, they can be seismically active for many 
reasons. Possible mechanisms for the generation of oscil- 
latio ns have bee n wide l y discussed i n the literature (see, 
e.g., iMcDermt^ il988h : IStergioul^ (l2003ll and references 
therein). Recently, much attention has been paid to r-modes 
- vortex oscillations that can be generated in rapidly rotat¬ 
ing neutron stars and accompanied by powerful gravitational 
radiation. In addition, oscillations can be excited in neu¬ 
tron stars, for example, during X-ray bursts (nuclear explo¬ 
sions in outermost layers of accreting neutron stars), burst¬ 
ing activity of magnetars (anom alous X-ray p ulsars and soft 
gamma-ray repeaters; see, e.g.. iKas^ 11200411 !. and glitches 
(sudden changes of spin periods) of ordinary pulsars. 

In this paper we focus of high-frequency (~ 100 kHz) 
pressure oscillations (p-modes) with high multipolarity I > 
100 localized in neutron star envelopes (cr usts). In our pre¬ 
vious paper iChugunov fc YakovlevlIioOFill we have studied 
these oscillations for I > 500. In that case p-modes are local¬ 
ized in the outer envelope (before the neutron drip point, at 
densities p < 4 x 10^^ g cm“®), where the equation of state 
(EOS) of stelar matter is relatively smooth. Accordingly, the 
oscillation spectrum is simple and well established. 

In the present paper we extend our analysis to p-modes 
with lower 1. These oscillations penetrate into the inner enve¬ 
lope of the star, where the EOS undergoes considerable soft¬ 
ening due to neutronization and becomes more complicated 
(essentially different for ground-state and accreted matter). 
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We show that the neutron drip affects strongly the oscilla¬ 
tion spectrum. If detected, this spectrum would give valu¬ 
able information on the EOS in neutron star envelopes and 
also on global parameters of neutron stars (their masses and 
radii). 


2 FORMALIZM 

Following IChueunov fc Yakovlev! we study oscilla¬ 

tions localized in a thin neutron star envelope. It is con¬ 
venient to use the approximation of a plane-parallel layer, 
and write space-time metric in the envelope as 

ds^ = dA — dz^ — B? (di?^ -I- sin^ i9dv5^), (1) 

where the local time t and local depth « are related to the 
Schwarzschild time t and circumferential radius r by 

t = t \Jl — Rg/R, z = {R —r)/\Jl — Ro/R, (2) 


r = R is the circumferential radius of the stellar surface, 
and (p are spherical angles, Rg ~ 2GM/c? is the gravitation 
radius, and M is the gravitational mass of the neutron star. 
The metric 0 is locally flat and allows us to use the New¬ 
tonian hydrodynamic equations for a thin envelope with the 
gravitational acceleration 


GM 

rVi - Rg/R' 


( 3 ) 


The pressure in the envelope is primarily determined by de¬ 
generate electrons and neutrons (in the inner envelope), be¬ 
ing almost independent of temperature T. Accordingly, we 
can use the same zero-temperature EOS for the equilibrium 
structure of the envelope and for perturbations. Employ¬ 
ing this EOS, we neglect the buoyancy forces and study p- 
modes. The linearized hydrodynamic equations (for a non- 
rot ating star) can be rewritten as (see, e.g., the monograph 

bylLa^Ei^ 

^ = c,^A<(. + <7.V<(., (4) 

where cj> is the velocity potential and Cg = dPo/dpo is the 
squared sound speed. The velocity potential can be pre¬ 
sented in the form 


cj, = e^‘^^Yi^{g,p)F{z), 


( 5 ) 


where lj is an oscillatio n frequency, and Yirr,. (tf, p] is a spher¬ 
ical function (see, e.g.. IVarshalovich et al.l lIlQSjjl l. An un¬ 
known function F{z) obeys the equation 


d^F 

d«2 c| da V c| 


R2 ) 


F = 0. 


( 6 ) 


The boundary condition at the stellar surface is F(Q) = 
0. It comes from the requirement of vanishing Lagrange vari¬ 
ation of the pressure at the surface. The formal condition 
lim^^oo F[z) = 0 in the stellar interior should be imposed 
to localize oscillations in the envelope. Of course, the actual 
variable a is finite and the real depth of oscillation localiza¬ 
tion will be controlled in calculations. The EOS of matter 
in neutron star envelopes contains a sequence of first-order 
phase transitions associated with changes of nuclides with 
growing density. These phase transitions are relatively weak 
(the density jumps do not exceed 20 per cent). We should 


add boundary conditions at all phase transitions within the 
envelope. These are two well known conditions at a plain 
boundary of two liquids The first condition 

can be written as 

F[{z) = FAz). (7) 

It ensures equal radial velocities at both sides of the bound¬ 
ary. The second condition is 

El = (8) 

pi \pi ) 

It comes from the requirement of pressure continuity at the 
boundary. Note, that the boundary conditions o and 
provide a source of buoyancy wh ich leads to the d ensity 
discontinuity of g-modes (see, e.g.. lMcDermottlll99fll . 

Oscillations of a plane-parallel layer for a polytropic 
EOS (P oc n being the p olytropic index) were stud¬ 

ied analytically bv IG oum (I199 Jl . In this case, the squared 
sound speed is Cg = pz/n. The solution for eigenfrequencies 
is 

+ A 

and eigenmodes are given by 

FAz) = exp [-^W+T) (2 ^W+V) I) , (10) 

where lI" (x) is a gen eralized Laguerre polynomial 
jAbramovitz fc StegunirL97 jf . and fc = 0, 1,... is the num¬ 
ber of radial nodes. 

Note, that the mode with k = 0 does not have any 
radial nodes; its properties are independent of the polytropic 
index n. This mode corresponds to the vanishing Lagrangian 
variation of the density (incompressible motion). Adding the 
condition A(j> = AU = 0 to Eq. 0, one can easily show that 
the mode with the frequency 

^ + 1) (11) 

and the eigenfunction Fo{z), defined by Eq. llOt . is the 
proper mode for a wide class of EOSs. Note, that the bound¬ 
ary conditions Q and are automatically satisfied for this 
mode, and it is continuous at phase transitions. The oscilla¬ 
tion frequency redshifted for a distant observer is 

= (^ “ ^) I 

= ^y/l-Ra/Ry^IifTT). (12) 

The frequency ujq will be used to normalize eigenfrequencies 
of other p-modes. The number of radial nodes k will be used 
to enumerate the modes. 


2.1 Self-similarity and scaling 

Let us use the equation of hydrostatic equilibrium dP/dz = 
pg and transform Eq. 0 taking the equilibrium pressure P 
as an independent variable. 


I 2 dF I ( l{l + A \r. n 

dP2 p (tp p2 g2 f,2 g2 p2 J 


( 13 ) 
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The boundary conditions 0 and 0 can be written as 


pi 


dFi 

dP 

Fi 


= p2 


dfa 
dP ’ 


- f2 +-II -^Pl -JW- 

pi \pi j dP 


(14) 

(15) 


Therefore, Eq. (O with the boundary conditions ITU , ira 
and with regularity requirement can be treated as the equa¬ 
tion for an eigennumber A = containing the scaling 

parameter = \/l{l + l)/(p-R) (with « l/{gR) for I 3> 1). 
Accordingly, the eigenfrequencies can be written as 


Wfc = g'^gkiO = Wo fkiO- 


(16) 


Here, pfc(C) and /fc(C) are functions which can be calculated 
numerically. They are universal for all neutron stars with a 
given EOS in the envelope. The velocity potentials Fk are 
also universal functions of P. Therefore, p-mode oscillations 
in stellar envelopes are self-similar and can be easily rescaled 
to a neutron star with any radius and mass. In principle, this 
can be used to determine R and M (see Sec. rra . 



Figure 1. The squared sound speed as a function of depth 
z. Solid line is for an accreted envelope. The dashed, dotted and 
dash-and-dotted lines are for the HP, RHS and SCM EOSa of the 
ground-state matter. 


3 NUMERICAL RESULTS 

Numerical results are presented for a “canonical” neutron 
star model, with the mass Me = 1.4M0 and the radius Rc = 
10 km. For this model, we have pc ~ 2.42 x 10^"* cm s“^, 

uJo « 1.56 X 10® [l{l + 1)/10'‘] s"^ (17) 

and (for a distant observer) 

wo = wo Vl - -Rg/P « 1.19x10® + s“\(18) 

Oscillation frequencies have been determined via a se¬ 
ries of iterative trials, checking the coincidence of the mode 
number and the number of radial nodes. 

3.1 Equations of state 

We employ two models of matter in neutron star envelopes, 
the accreted and grou nd-state matter. For the a ccreted mat¬ 
ter, we use the EOS ol lHaensel fc ZdunikI ilOPOll (HZ). It was 
derived by following transformations of atomic nuclei (beta 
captures, emission and absorption of neutrons, pycnonuclear 
reactions) in an accreted matter element with increasing the 
pressure. The EOS was calculated for the densities from 
p — 3.207X 10^ g cm~® to 1.462x10’^® g cm“®. For lower den¬ 
sities, we have taken the matter composed of ®®Fe and the 
EOS of degenerate electrons with electrostatic corrections. 
For higher densities, we use the EOS of the ground-state 
matter presented by iBavm. Pethick fc S utherland) ^W7^ 
(BPS) because, as remarked bv iHaensel fc ZdunilJ iImQcF . 
the HZ EOS becomes very similar to the BPS EOS at 
p > 10®® g cm“®. 

We have also considered envelopes composed of the 
ground-state (cold catalyz ed) matter. In th e outer enve¬ 
lope we use the EOS of |HaenseL£Mhcho3 ll994| ) (HP) 
and th e recent EOS of lP.iister. Hemnel fc Schaffner-Bielichi 
l2QQ5|) (RHS). For the inner envelope, we employ the EOS 
of iNegele fc VautherirJ (Il973ll . 

Phase transitions in these EOSs have been treated care¬ 
fully using the boundary conditions 0 and 0 at any phase 


transition. For comparison, we have also employed the model 
of the outer envelope composed of ground-state matter with 
a smoothed composition (the smooth composition model - 
SCM). In the latter case we have included only a large den¬ 
sity jump at the neutron drip boundary between the inner 
and outer envelopes. 

The squared sound speed cf as a function of depth 2 
for all these EOSs is shown in Fig. Q The solid line is for 
the accreted envelope; the dashed, dotted and dash-and-dot 
lines are for the HP, RHS and SCM EOSs of the ground-state 
matter. The different versions of the ground-state EOS show 
approximately the same sound speed profiles, but the profile 
in the accreted envelope is significantly different. 

The depth of the accreted envelope (up to the density 
2.004 X 10®"® g cm“®, which is the largest density in the 
envelope, where the atomic nuclei are present, for the BPS 
EOS) is 2 « 1150 m. For all models of the ground-state 
matter, the largest density in the envelope has been taken 
« 1.7 X 10®"® g cm“®; the envelope depth is 2 « 985 m. 

3.2 Eigenfrequencies 

Figures]^ and 13 show squares of dimensionless eigenfrequen¬ 
cies uf/uiQ versus multipolarity I for accreted and ground- 
state envelopes of the canonical neutron star. Because of the 
scaling l|T^ the figures can be easily transformed to a star 
with any gravity g and radius R by changing scale of the I 
axis by a factor of g R/[g^ Re)- 

For any envelope, the modes with I > 300 can be subdi¬ 
vided into two groups, with a pronounced linear dependence 
and with a weak dependence of on Z. As will be shown 

in Section f3.4l the modes of the first type [the inner modes, 
shown by thicker lines in Fig. 0 are localized in the vicinity 
of the neutron drip point, while the modes of the second 
type {the outer modes) are localized in the outer envelope. 
In Figs. 0 and 0 one can see a number of quasi-intersections. 
When passing through a quasi-intersection point (with grow¬ 
ing 1), an inner mode gains an additional radial node, but 
an outer mode loses one. 


















4 A. I. Chugunov 


CM O 

3 


6 

5 
4 

3 
2 
1 

0 

100 200 300 400 500 600 700 800 900 1000 

I 


ri ■ ■ 1 1 ■ 1 1 ■ 1 1 1 1 1 1 1 1 ■ 1 1 ■ 1 1 ■ 1 1 ■ 1 1 1 1 1 1 1 1 ■ ■ 1 ■ 1 1 1 ■ 1 

6 

ri ■ ■ 1 1 ■ 1 1 ■ 1 1 1 1 1 1 1 1 ■ 1 1 ■ 1 1 ■ 1 1 ■ 1 1 1 1 1 1 1 1 ■ ■ 1 ■ 1 1 1 ■ 1 j 

■ Accreted envelope ® . 


■ Model envelope 6 .; 

^ 5 . / 

5 

: . ^ 1 

/■.... .. 4 ...U 

4 

^ .. ../ 

..T. 3 

CM 0 

3 

CO 

\ 

\ 

\ 

i 

1 


0 

CM ^ 

3 

-- ■ 2 

- 1 - 

2 

- . 1 

0 - 

1 

0 

’ 0 


: 1 


100 200 300 400 500 600 700 800 900 1000 

I 


Figure 2. Squared normalized eigenfrequencies versus 

multipolarity I for the accreted envelope of the canonical neutron 
star. The numbers next to curves indicate the number of radial 
nodes. Thin parts of the curves correspond to the outer modes 
and (for low 1) modes which a spread over the entire envelope, 
while thick segments refer to the inner modes. 
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Figure 3. Same as in Fig.l2| but for the envelope composed of the 
ground-state matter (the inner modes are not emphasized). Lines 
are for the HP EOS; crosses ‘x’ are for the RHS EOS; crosses ‘H-’ 
are for the SCM. 


Let us consider the outer modes. The eigenfrequen¬ 
cies are the same (within ~ 1%) for all ground-state EOSs 
(see Fig. 0 . For the accreted envelope, eigenfrequencies are 
larger because the EOS is stiffer. With decreasing I, oscil¬ 
lations penetrate deeper into the outer envelope, where the 
EOS is softer because electrons become relativistic, and be¬ 
cause they undergo beta-captures. It leads to a gradual de¬ 
crease of loIIloq. As in the model with the polytropic EOS, 
given by Eq. 0 , separations between squares of neighbor¬ 
ing dimensionless eigenfrequencies ujI IlJq are approximately 
constant for a fixed 1. The weak decrease of separations with 
the growth of k is due to the penetration of oscillations into 
deeper layers of the star, where the EOS is softer. The latter 
effect is more pronounced for the ground-state matter ow¬ 
ing to a stronger softening of the EOS. Einally, at I ~ 1000 
the outer p-modes are localized in the outer layers of the 


Figure 4. Same as in Fig.|21 but for a model envelope in which 
the outer envelope is composed of the ground-state matter and 
the inner envelope is composed of ^'^®Se. 


outer envelope, where the matter is composed of ®®Fe nu¬ 
clei for both accreted and ground-state EOSs. Accordingly, 
eigenfrequencies become nearly equal. 

Naturally, the oscillation frequencies of outer modes 
with I > 500 f or the ground state envelope are the same 
as calculated bv lOhiigunov fc Ya.kovle^ || 22 Q^. 

To demonstrate explicitly that inner modes are caused 
by the neutron drip in the inner envelope, in Fig. |1| we 
present eigenfrequencies for a model envelope without any 
neutron drip. Here we employ the RHS EOS in the outer 
envelope but assume that the inner envelope is composed of 
^^®Se ions (the last element at the outer envelope) and elec¬ 
tron gas (no free neutrons). The oscillation spectrum does 
not contain any inner modes. A small decrease of loIIloq for 
200 < I < 500 is produced by the softening of the EOS at 
the bottom of the outer envelope. The growth of frequen¬ 
cies at I ~ 100 is caused by the penetration of oscillations 
into the inner envelope, where our model EOS is polytropic 
(with the index n = 3). Accordingly, oscillation frequencies 
tend to the values provided by the polytropic model 0 . 


3.3 Inferring M, R, and the crustal EOS from 
oscillation spectrum 

If detected, outer modes would give us uio, and there¬ 
fore M R~^ ~ Rg/R', see Eq. I12II . A detection of the 

only one fundamental mode {k = 0 ) would be sufficient 
to determine M R~^ -^1 — Rg/R- A detection of several 
outer modes (with different I and/or k) would confirm and 
strengthen this determination. 

Our calculations show that for the inner modes the ratio 
wfn/tJo is a linear function of 1. Using the scaling relation 
we can present this linear dependence in the form 

~ ^ B I, B = P/{gi4Re)- (19) 

Here, <714 is the surface gravity in units 10^"* cm s“^. Re — 
77/10® cm=77/10 km, while A and /3 are dimensionless con¬ 
stants determined by the EOS in a neutron star envelope. 

Eor the canonical neutron star with the ground-state 
envelope, we obtain A — 0.75 and B = 0.0032 in the case of 
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Figure 5. Angle-averaged energy density of oscillations for 
modes with I = 100 in the accreted envelope. The subscript of 
£ indicates the number of radial nodes. The root-mean-square 
amplitude of radial displacements at the stellar surface is 1 m. 
The vertical dotted line shows the boundary between the inner 
and outer envelopes. 

inner modes with lowest freqnencies. For the same star with 
the accreted crnst we have A = 0.65 and B = 0.0073. The 
values of B allow us to determine /3. In this way we obtain 

A = 0.75, /3 = 0.0013 for gronnd state crnst; (20) 

A = 0.65, /3 = 0.0030 for accreted crust. (21) 

Hence, the difference between the ground-state and accre¬ 
tion envelopes is quite pronounced in oscillation spectra. 

Therefore, if several (minimum two) inner modes could 
be detected in addition to outer modes, their frequencies 
could be htted by a function and the values of A and 
B could be determined. An accurate determination of A 
would enable one to distinguish between the ground-state 
and accretion envelopes. The value of B would give then gR. 
Combining this gR with the value g R~^ ^1 — Rg/R, deter¬ 
mined from the detection of the outer modes, one would get 
a simple system of two equations for two unknowns, M and 
R. Thus, a detection of one outer mode and several inner 
ones could in principle enable one to discriminate between 
the ground-state and accreted envelopes and determine neu¬ 
tron star mass and radius. 


3.4 Eigenmodes 

Figures 01^ show profiles of the angle-averaged energy den¬ 
sity of oscillations as a function of 2 . The root-mean-square 
amplitude of radial displacements of the stellar surface has 
been set equal to 1 m. The subscript of e indicates the 
number of radial nodes. The vertical dotted line marks the 
boundary between the inner and outer envelopes {z ~ 432 m 
for the accreted envelope, and 2 « 364 m for all ground-state 
envelopes of the canonical neutron star). 

Figure depicts eigenmodes with I = 100 for the ac¬ 
creted envelope. The modes are spread over the entire enve¬ 
lope; their subdivision into the outer and inner modes is not 
obvious. However, some traces of two mode types are visi¬ 
ble. The mode with one radial node (the dashed line), whose 



Figure 6. Same as in Fig. for modes with I = 400 in the 
accreted envelope. 


frequency belongs to the branch of the outer modes, is pri¬ 
marily localized in the vicinity of the neutron drip point. 
However, other modes do not demonstrate this feature. The 
effects of phase transitions are relatively small (~ 20%) and 
only slightly noticeable in Figure]^ They are local and do 
not change global (on scales of > 20 m) energy density pro¬ 
files. 

Figure El shows eigenmodes with I = 400 for the ac¬ 
creted envelope. This value of I is very close to the quasi¬ 
intersection point for the modes with 4 and 5 radial nodes 
(see Fig. EJ. The subdivision into the outer and inner modes 
is clear - the modes with k =0, 1, 2, 3, 5, 6 radial nodes 
are localized in the outer envelope, but the energy of the 
mode with k =4 is concentrated in the outer part of the in¬ 
ner envelope. This subdivision is the same as in Section f3.2l 
(see FigureEJ. The energy-density profiles of the fourth and 
fifth modes are very similar at 2 < 250 m, but at 2 ~ 500 m 
the energy densities differ by more than two orders of mag¬ 
nitude! The outer modes ’feel’ the lowering of the sound 
speed in the outer layers of the inner envelope (see Fig. 0 
and increase their energy density in this region. However, 
the increase is not so large as for the inner modes. The sig¬ 
natures of phase transitions are very small (^ 10%) and 
hardly visible in Figure El They are local and do not change 
energy density profiles on scales > 20 m. 

Figures El and El are plotted for the ground-state enve¬ 
lope with the RHS EOS. The results for the HP and SCM 
are qualitatively the same. 

Figure El depicts eigenmodes with I = 100. The modes 
are localized in the entire envelope and cannot be subdivided 
into the outer and inner ones. The traces of these two types 
of modes are weaker than for the accreted envelope (see Fig. 
EJ. The signatures of phase transitions in the outer envelope 
are small < 10% and have scales ~ 10 m. They are noticeable 
only for the modes with a few number of radial nodes. Many 
modes show large (~ 50%) jumps of the energy density at 
the neutron drip point. 

FigureElshows eigenmodes with I = 500. This value of I 
is close to the quasi-intersection point for modes with k =2 
and 3 and with k =5 and 6 (see Fig. EJ- The modes can 
obviously be subdivided to the two types. The modes with 
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Figure 7. Same as in Fig. [3 for modes with I = 100 in the 
ground-state envelope. 



Figure 8. Same as in Fig. for modes with I = 500 in the 
ground-state envelope. 


k =0, 1, 3, 4, 5 are localized in the outer envelope; the en¬ 
ergy of the modes with k =2 and 6 is concentrated near the 
neutron drip point. The subdivision of modes is the same as 
suggested in Section Em on the basis of Figure El The en¬ 
ergy profiles of the second and third modes are very close for 
2: < 200 m, but for z ~ 420 m the energy density differs by 
more than three orders of magnitude. Qualitatively the same 
feature is demonstrated by the fifth and sixth modes. The 
outer modes ‘respond’ to the lowering of the sound speed af¬ 
ter the neutron drip point by increasing the energy density 
in this region. This increase is greater for the second and 
fifth modes whose frequencies are close to the frequencies of 
the inner modes. However, it is not so large as for the in¬ 
ner modes. The signatures of phase transitions in the outer 
envelope are extremely small 5%) and are almost invisi¬ 
ble in Figure El The scales of such features are ~ 1 m. The 
large phase transition at the neutron drip point produces 
the signature with the same properties. 


4 CONCLUSIONS 

We have studied high-frequency pressure oscillations which 
are localized in the envelopes of neutron stars composed of 
the accreted or ground-state matter. 

Our main conclusions are as follows. 

(1) The oscillations are almost insensitive to various 
modifications of the EOS for the ground-state matter (sec¬ 
tion EHJ. All EOSs we have used (HP, RHS, SCM) give the 
same oscillation spectrum. 

(2) The neutron drip and associated softening of the 
EOS in the inner envelope do not affect strongly the spec¬ 
trum of the well known (outer) oscillation modes which are 
localized predominantly in the outer envelope (section EHJ. 

(3) However, the neutron drip leads to the appearance 
of inner oscillation modes localized mostly near the neutron 
drip point (sectionEJ. The spectrum of these modes is sensi¬ 
tive to the EOS in the envelope (accreted or ground-state). 

(4) The p-mode oscillation problem is self-similar (in the 
plane-parallel approximation). Once the problem is solved 
for one stellar model, it can easily be rescaled to neutron 
star models with any mass and radius (but the same EOS 
in the envelope; see section lOl . 

(5) A detection and identification of one outer mode 
and several inner modes would enable one to discriminate 
between the ground-state and accreted envelope and deter¬ 
mine neutron star mass and radius lsection l3.3t . For exam¬ 
ple, a detection of the fundamental mode with I = 900 at the 
frequency 74 kHz and of two inner modes with I — 300 and 
I = 900 at 56 kHz and 140 kHz, respectively, would indicate 
a canonical neutron star with the ground-state envelope. 

Therefore, high-frequency pressure modes are excel¬ 
lent tools to explore the physics of matter in neutron star 
envelopes and to determine masses and radii of neutron 
stars. The oscillation frequencies could be detected by radio- 
astronomical methods very precisely. A detailed analysis 
of pulse shapes of some radio pulsars reveals that oscil- 
lations with large multi polarity are possibly excited there 
iClemens fc R,osenll2004^ but their frequencies are ~30 Hz, 
so that they are not high-frequency p-modes we discuss here. 

A search for high-frequency p-modes could be useful. 
High-multipolarity p-modes do not damp very quickly be¬ 
cause they do not produce any po werful gravitational or 
electr omagnetic emission (see, e.g., IChugunov fc Yakovlev! 
[200^. They are robust because they are relatively indepen¬ 
dent of the thermal state of the star, and they should not 
be strongly affected by neutron star magnetic fields. The in¬ 
ner p-modes, localized in the inner envelope, could be easily 
triggered by pulsar glitches, which are thought to occur just 
in inner envelopes of pulsars. IChugunov fc Yakovlev! (l2005h 
studied the dissipation of p-modes localized in the outer 
envelope; this dissipation is mainly produced by the shear 
viscosity. It may be enhanced by thin viscous layers near 
numerous nuclear phase transitions. The viscosity in these 
layers can be diffusive or turbulent. Note, that fundamen¬ 
tal modes do not produce viscosity layers because they pass 
phase transitions without velocity discontinuities (see Sec. 
EJ Their dissipation is not enhanced by phase transitions. 

Finally, p-modes in neutron star envelopes are rela¬ 
tively insensitive to the EOS and composition of neutron 
star cores. However, these modes can be useful to discrimi¬ 
nate between ordinary neutron stars and strange stars with 
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crust. The latter stars are thought to contain extended cores 
composed of strange quark matter. Nevertheless, a core is 
assumed to be surrounded by an envelope of normal mat¬ 
ter (see, e.g., |zd uni ^l2nn2h . so that a strange star with the 
crust may look like an ordinary neutron star from outside. 
The density of the normal matter in a strange star does not 
exceed the neutron drip density. The pressure modes in the 
envelopes of such stars should easily ‘feel’ underlying dense 
quark matter, and the oscillation spectrum would reflect the 
presence of the quark core. We intend to consider this effect 
in a future publication. 
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